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Abstract 

We have investigated the electronic structure and magnetic properties of GdN as a function 
of unit cell volume. Based on the first-principles calculations of GdN, we observe that there is 
a transformation in conduction properties associated with the volume increase: first from half- 
metallic to semi-metallic, then ultimately to semiconducting. We show that applying stress can 
alter the carrier concentration as well as mobility of the holes and electrons in the majority spin 
channel. In addition, we found that the exchange parameters depend strongly on lattice constant, 
thus the Curie temperature of this system can be enhanced by applying stress or doping impurities. 

PACS numbers: 75.10.-b, 71.27.+a, 71.70.Gm 
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Electronic and transport properties of the rare-earth nitrides have long been a challenge 
to investigators: the nitrides are difficult to fabricate into single phase crystals and the 
experimental picture of their electronic structures is far from clear. Although most rare- 
earth nitrides have shown to be semi-metallic, there are still uncertainties about GdN. Based 
o. the d.ect ..es.t,vitv ,„ea..e.ents. X,ao . al. concluded that GdN . an insulate. R. 
Nevertheless, the values they obtained for the resistivity do not rule out the possibility of 
a semi-metallic state. Furthermore, there is no clear-cut picture emerging from a series of 
stud,es by Kald,s and co-autho. in 70. and SO. flBQ. They pointed out that the specific 
resistivity decreased with increasing temperature, which suggests that the GdN system is 
semiconducting P], while their later experiments on better quality samples showed a very 
arge carrier concentration of 1.9 x 10^^ cm~'^, which is more characteristic of a semi-metal 
^. Even though the optical absorption edge for GdN is about 1 eV ^, which is similar 
to other rare-earth nitrides, at issue is the position of the chemical potential relative to the 
band edges in both the spin majority and spin minority band structures. 

An appealing property of GdN is that it is ferromagnetic with a large gap at the Fermi 
energy in the minority spin states, according to the electronic structure calculations based 
on the local density approximation (IDA) QSQ. At the same time GdN is semi-metallic 
in majority spin states with electron and hole pockets at the Fermi surface |6j. This latter 
propertyhas led to some interest in GdN as a possible candidate for spin-dependent transport 
devices js| , exploiting the spin filter, giant magnetoresistance or tunneling magnetoresistance 
effects. 

An accurate description of the electronic structure of rare-earth compounds is a very 
challenging problem because of their unffiled 4/ shells jol. Calculations based on local spin 
density approximation (LSDA) are well known to underestimate the band gap in semicon- 
ductors. Thus LSDA and similar computational methods may not be able to correctly 
describe whether a highly correlated system, like GdN, is semi-metallic or semiconducting 
at the cuiUhHun, volume a. Nonetheless, ,f we are interested in the trend of how the elec- 
tronic and magnetic properties vary with the change of volume, we can obtain a reasonable 
picture for GdN from LSDA with additional Hubbard correlation terms describing on-site 
electron-electron repulsion associated with the 4/ narrow bands (LSDA+f/ approach) jlol |. 
Actually, due to the fact that the / states of this system are exactly half occupied, there is 
no orbital moment and the anisotropic and multipole effects are minimal. As a result, GdN 
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is the ideal material to study the magnetic exchange interactions in rare-earth nitrides ^ | . 

In this letter, we show that applying stress can influence signiflcantly the electronic and 
magnetic properties of GdN. Using the flrst-principles approaches, we demonstrate that the 
system exhibits nominal "half-metallic" band structure at the equilibrium lattice constant, 
and then semi-metallic and/or semiconducting character develops with increasing lattice 
constant. We note that the magnetic properties are also extremely sensitive to the volume 
variations, i.e., the exchange interactions are at flrst ferromagnetic, then the calculated 
magnitudes of exchange parameters reduce substantially with increasing volume, suggesting 
that the Curie temperature is reduced with an increase in lattice constant. 

The flrst-principles band structure approach applied in this work is the full-potential 



12| . In the total energy calculations. 



linear-augmented-plane-wave plus local-orbital method 
the factor RMrKmax is chosen to be 8. We found that using as large as 4000 fc-points in 
the Brillouin zone was necessary to obtain the energy convergence up to 0.1 meV. Following 
the previous work [3], we used a Hubbard U = 6.7 eV and an exchange J = 0.7 eV for Gd 
compounds in the LSDA+?/ scheme The calculated density of states (DOS) of GdN 
agrees well with photoemission data of nitrogen covered Gd(OOOl) surfaces as shown 
in Fig. H Given that both experiment and theory share the similar lattice structures, this 
suggests that our computational approach correctly represents the electronic structure of 
GdN. 

When we studied the band structure, we found that, for the theoretical lattice constant 
a = 4.92 A [Fig. Efa)], our LSDA+f/ calculations show that there exist a hole pocket at F 
point and an electron pocket at X point. Strong hybridization between Gd 5d and N 2p spin 
majority states is clearly shown around the X point, while there is no such hybridization 
for spin minority states, indicating GdN to be half-metallic, with a gap about 0.6 eV in the 
spin minority channel. As we increase the lattice constant, this system gradually develops 
semi-metallic features in the majority states: e.g., the d — p hybridization around X point 
disappears and the Fermi level down shifts to the top of the hole section at F point [Fig. 
Efb)]. With further increase of volume the system eventually becomes semiconducting [Fig. 
|2fc)]. The exact volume at which metal-semimetal or semimetal-semiconductor transitions 
occurs might not be precisely determined from LSDA-based calculations. However, the trend 
is clear and the predicted transitions could be observed experimentally. 

We note that volume expansion or applied strain can be used to control the carrier density 
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as well as to some degree the mobility of the carriers, as can be seen from the band dispersion 
curves near the Fermi energy, in Fig. [21 At larger volumes the electron (hole) pockets become 
substantially shallower, meaning that the bottom (top) of the band approaches the Fermi 
energy and the area of the Fermi surface decreases. This is quite clear, for example, when 
comparing the hole pockets around F point in Fig. |2ta) and (b), especially those Fermi- level 
crossing points. As a result, the density of states at the Fermi energy reduces by about 50% 
when lattice parameter is increased by only 5%. And the decrease of dEk/dk (tangent to the 
band dispersion) at Fermi level can also be seen from Fig. indicating that the mobility of 
the holes decreases with the increase of the cell volume. 

Furthermore, as also can be seen from Fig. |2l the indirect band gap of GdN can be 
modified by hydrostatic pressure. Based on the prediction of electron-hole-liquid theory 
[iS^ . a first-order semiconductor to semimetal transition starts to take place as the indirect 
band gap decreases with decreasing lattice constant. Thus it is conceivable that with proper 
control of the indirect gap via external pressure, it is possible to explore the ground state of 
this correlated electron-hole liquid, and tune experimentally the semimetal phase transition. 

We have also explored the modifications of the band structure and electronic properties 
due to the biaxial strain, which can be produced by epitaxially growing the film on the 
substrate with a different lattice constant. When the lattice parameter of the substrate is 
different from that of equilibrium GdN, the lattice mismatch can produce strain. In such 
cases, the lattice symmetry of GdN becomes tetragonal with c/a ratio less than 1 (larger 
substrate lattice constant, tensile strain) or larger than 1 (smaller substrate lattice constant, 
compressive strain). We performed total energy calculation on the Poisson ratio u for GdN 
by applying the biaxial strain in ab plane and observing the variation of c, then by comparing 
the biaxial strain given as 

Ac/ Co 2v 
Aa/aQ 1 — z/' 

we deduced a z/ = 0.2, which is small compared to metals but similar to the known value of 
TiN. We found that the trend of band structure change, due to the biaxial strain, remains 
qualitatively the same as that caused by volume strain, which is expected because of the small 
Poisson ratio. The unit cell volume of GdN increases due to biaxial tensile stress, hence the 
modification of the band structure near the Fermi energy is similar to that observed during 
volume expansion. This may imply that the tensile stress tends to force the system to be 
less metallic and decreases the density of states near Fermi surface. When applying the 
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compressive strain, the system tends to become more metallic. Based on our calculations, 
the minority spin shows considerable reduction in the band gap with bands cross over Fermi 
energy with in-plane compressive strain of 3% and then the system is no longer half-metallic. 

We also found that the change in electronic structure with the applied stress significantly 
affects the magnetic properties of GdN, particularly, the exchange interactions. We analyzed 
the exchange interactions using Heisenberg Hamiltonian: 

H = J n Si ■ Sj , 

n i>j 

where J„ are exchange parameters and n is the index of the nearest-neighbor shell. In these 
calculations, we limited our considerations to the third nearest neighbor interactions, i.e., n 
runs from 1 to 3. The exchange parameters used in the model Hamiltonian, are obtained 
from the first-principles band structure calculations. To do this, we carried out total energy 
calculations on four different magnetic ordering configurations of the fee structure. One is 
ferromagnetic (FM) ordering, and the other three are anti-ferromagnetic (AFM) orderings, 
as have been described in Ref. and here for convenience are called AFM I, II and III, 
respectively. According to the model Hamiltonian, the total exchange energies per magnetic 
lattice site of the four magnetic orderings can be explicitly expressed as: 

Efm = Eo + 6Ji + 3J2 + 12 J3, 
Eafmi = Eq — 2Ji + 3J2 — 4J3, 
EafMii = Eq — 3J2, 

EaFMjii = Eq — 2Ji + J2 + 4J3, 

where Eq is the reference energy. Based on these equations and the results obtained from the 
total energy calculations, those J parameters can be deduced accordingly. Due to the fact 
that the energy differences between these ordering states are generally very small, extreme 
care is needed in the calculations. Hence, comparisons are only made between FM and AFM 
energies calculated on the same structure and with the same computational parameters, to 
avoid any error caused by the different symmetries or shapes. 

Our LSDA+ U calculation gives the correct ground state for GdN, i.e., FM ordering. The 
Jn values are, nonetheless, quite small, which is expected from the low transition tempera- 
ture. In addition, we found that these exchange parameters of GdN depend strongly on the 
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lattice constant (Fig. E)). When the lattice constant is increasing, the FM Ji decreases and 
J2 even changes sign, namely from FM to AFM, while J3 keeps almost the same. Thus the 
FM transition temperature T^, which is proportional to the sum of neighboring exchange 
energies 12Ji + 6J2 + 24J3 according the mean-field theory, is sensitive to the change of 
lattice constants. 

It is well known that the oscillatory Rudermann-Kittel-Kasuya-Yosida-type interaction 
varies sensitively with the density of charge carriers. Thus, based on the above observations, 
we can see that the enhancement of the exchange interactions between neighboring magnetic 
sites, when the lattice constant decreases, is caused by the increase of the number of free 
charge carriers as the GdN system becomes more metallic. Actually, the strong lattice con- 
stants dependence of the J parameters of GdN is a manifestation of strong lattice constants 
dependence of its electronic structure, which lies between the metal and insulator phase. 
Furthermore, the trend of J2 with the increase of the volume implies the strengthening of 
an AFM superexchange interaction when the system becomes less metallic. This can be 
understood as a competition between the superexchange and indirect carrier mediated ex- 
change interactions. Usually the former is roughly proportional to t'^/U, where t is the band 
energy or hopping integral, and does not change too strongly when the lattice parameter 
increases. Whereas the later, as we already discussed, would decrease with the increase of 
the lattice constant. Hence this reduction causes an overall increase in the antiferromagnetic 
coupling between second-nearest-neighbor Gd sites. Therefore we can see that the magnetic 
properties of Gd nitride are strongly related to their electronic properties. Based on these 
findings, we expect that GdN, in which a strong AFM and FM competition exists, could 
be AFM ^ or even more complicated structure such as spin-glass [17J when experiencing 
different stress. 

Monte Carlo (MC) simulations, based on the model Heisenberg Hamiltonian with ab initio 
derived exchange parameters are used to obtain the Curie temperature. Same method has 
been applied successfully in the study of complex permanent magnetic materials The 
lattice studied in our MC simulation is a 10a x 10a x 10a fee cell (4000 spins) with periodic 
boundary conditions, where a is the lattice constant. At theoretical lattice parameter (a = 
4.92 A) we find the Curie temperature to be about 38 K, agreed reasonably well with the 
experimentally observed 58 K. This agreement is quite impressive, considering the strong 
correlated nature of this system. One possible reason for the underestimate of Tc may arise 
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from neglecting the correlated hopping processes in one-electron picture, which are higher 
order processes proposed by Kasuya and Li to explain of ferromagnetic exchange in GdN 
[19I 1 . In addition, we found that the energy levels of 4/ states are crucial in determining the 
accurate exchange parameters. Our calculation on equilibrium structure of GdN, based on 
the present U parameter, gives the unoccupied 4/ states 5 eV above Ep and occupied 4 / 
states 6 eV below, which is in good qualitative agreement with the experimental situation 
j^l, providing support for the validity of our studies of the magnetic properties of GdN. 

Ordinarily, we would expect the ferromagnetic ground state to provide a stabilizing role 
to the spin- dependent electric conduction. Thus to enhance the Curie temperature of GdN 
to practical range is quite important. As we already mentioned above, one possible way 
is to apply stress. Moreover, doping impurities can also be an alternative: it was reported 
that by adding few percent of magnetic ions such as Mn, the Tc of GaN or ZnO well exceeds 



room temperature range 



21[, hence we believe this is a promising remedy for concern that 



the Curie temperature of GdN is on the low side. 

In summary, we have found that there is a large lattice constant dependence of the elec- 
tronic and magnetic properties in the spin dependent band structure of GdN. This material, 
with half-metallic gap about 0.6 eV, exhibits a ferromagnetic ground state, rendering it 
an attractive candidate for spintronic devices. Hence we appeal strongly for experimental 
efforts to study this interesting compound and validate our claim. 
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FIG. 1: Comparison between the calculated DOS (majority spin: solid line, minority spin: dotted 
line) at theoretical lattice constant (a = 4.92 A) and the photoemission spectra for nitrogen on 
Gd(OOOl). The new photoemission features that are not attributable to the Gd(OOOl) substrate 
compare well with the calculated DOS for GdN. Features with strong N 2p weight are indicated. 
Binding energies for experiment are shifted to higher binding energies as expected with a final state 
spectroscopy of a correlated electron system, but the shift is roughly uniform for key features of 
the photoemission spectra (taken for both s and p polarized light). 

FIG. 2: Band structure of GdN in the vicinity of the Fermi energy for three volumes: (a) at the 
calculated equilibrium lattice parameter a = 4.92 A, (b) at the lattice parameter increased by 5% 
(a = 5.16 A), (c) at the lattice parameter increased by 14% (a = 5.63 A). Solid and dotted lines 
represent spin majority and spin minority states, respectively. The change of conducting properties 
are indicated by the change of energy difference between the top (bottom) of the hole (electron) 
pockets and the Fermi energy. 

FIG. 3: Exchange parameters of GdN plotted as function of the lattice constant; notice the strong 
volume dependence of Ji and J2. 
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